function [ coeficiente_determinacion ] = CoefDeterminacion( x,y )% x es la recta pendiente y es el balance
% x es la recta pendiente y es el balance

% p = polyfit(x,y,1);
% yfit = polyval(p,x);
yfit=x;
yresid = y - yfit;
SSresid = sum(yresid.^2);
SStotal = (length(y)-1) * var(y);
coeficiente_determinacion = 1 - SSresid/SStotal;
%coeficiente_determinacion = 1 - SSresid/SStotal * (length(y)-1)/(length(y)-length(p)-1);
 %coeficiente_determinacion=corr2(x,y);
 
%  c=corrcoef([x,y]);
%  r=c(1,2);
%  coeficiente_determinacion=r^2;

end

